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Abstract 

A preliminary measurement of the branching fraction B(t~ — » Kgir~u T ) is made using 384.6 fb -1 
of e + e~ collision data provided by the PEP-II collider, operating primarily at i/i = 10.58 GeV, 
and recorded using the BABAR detector. From this we measure: B(r~ — ► K Qi k~v t ) = (0.840 ± 
0.004 (stat) ± 0.023 (syst))%. This result is the most precise measurement to date and is consistent 
with the world average. 
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1 INTRODUCTION 



Hadronic r decays provide a clean laboratory for studying the hadronic weak current. Hadronic 
products from r decays give access to the light quark vector (V) and axial-vector (A) spectral 
functions, which give insight into the dynamics of QCD at intermediate scales as well as providing 
tests of the Standard Model itself pQ. 

For hadronic r-decays, SU (3) t symmetry breaking can be used to determine the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element magnitude \V US \ J2], the strong coupling constant, a s , 
and the strange quark mass, m s [3]. These are all tests of the Standard Model as deviations from 
values measured in other processes would indicate new physics. 

Hadrons from r decays are produced via W emission. Relative to non-strange (ud) currents, 
strange (us) currents of r decays are suppressed by (|V us |/|T4d|) 2 — tan 2 # c , where \V u d\ and \V US \ 
are the absolute values of the CKM matrix elements. Resonant decay dominates these currents: 
the strange vector current is dominated by a K* resonance which decays to Kir and the strange 
axial-vector current by the K\ which decays mostly via Kp and K*n to Knir. 

The high luminosity provided by the PEP-II collider, coupled with a large t + t~ production 
cross-section near the operating energy of = 10.58 GeV, provides a large data sample with which 
to study the strange hadronic decay r~ — > (Ki:)~v T using the BABAR detector. Measurements of 
r decay branching fractions to strange hadronic final states and studies of the strange spectral 
functions have been conducted by ALEPH [3], CLEO [5] and OPAL [6], but have been limited by 
statistics. About a hundred times larger samples of r-events have been provided by the 5-factories 
BABAR and BELLE. The measurement of B(t~ — > KgTT~v T ) was recently carried out by BELLE 
[7] with a significantly reduced uncertainty compared to that of previous measurements. 

2 THE BABAR DETECTOR AND DATASET 

The BABAR detector is described in detail in [8] . Charged particles are detected and their momenta 
measured with a 5-layer double sided silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH) inside a 1.5 T superconducting solenoidal magnet. A ring-imaging Cherenkov detector 
(DIRC) is used for the identification of charged particles. Energies of neutral particles are measured 
by an electromagnetic calorimeter (EMC) composed of 6, 580 CsI(Tl) crystals, and the instrumented 
magnetic flux return (IFR) is used to identify muons. 

The analysis described in this paper is based on data taken using the BABAR detector at the 
PEP-II collider [9] located at SLAC in the data-taking periods between October 1999 and August 
2006. During this period a total of 384.6 fb _1 of data was recorded with a cross-section for t + t~ 
pair production of (0.919 ± 0.003) nb [15]. This data sample contains over 700 million r decays. 
Monte Carlo (MC) studies of simulated signal and background events were carried out using various 
MC samples. The r MC events studied were generated with KK2f [10J and decayed with TAU0LA [11] 
using r branching fractions based on PDG 2006 [12J. In the MC, the r~ decays to Kgir~ via the 
K*(892)~ resonance with a branching fraction of 0.90%. Non-r hadronic and dilepton MC samples 
are used for studying the non-r backgrounds. 

3 SELECTING r K%-K~v r EVENTS 

In this analysis events containing a pion and a Kg in the final state are studied, where the K s is 
reconstructed in the 7r + 7r~ mode. 
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The event is divided into two hemispheres in the center-of-momentum system (CMS) using the 
plane perpendicular to the thrust axis, which is the direction which maximizes the sum of the 
longitudinal momenta of the neutrals and tracks in the event [16j . One hemisphere of the event 
is required to contain only one charged track; this is defined as the tag hemisphere. The other 
hemisphere is required to contain three charged tracks; this is defined as the signal hemisphere. 
The tag track and at least one of the signal hemisphere tracks are required to point towards the 
interaction point. 

Approximately 35 % of r decay to fully leptonic final states. Requiring the track in the tag 
hemisphere to be identified as an electron or muon while requiring the signal hemisphere to contain 
only hadrons strongly reduces backgrounds from e + e~ — ► qq events. Electrons are identified using 
the ratio of calorimeter energy to the track momentum (E/p), the ionisation loss in the tracking 
system (dE/dx) and the shape of the shower in the calorimeter. Muons are identified by hits in 
the IFR and small energy deposits in the calorimeter. 

Kg candidates are constructed from any two oppositely charged tracks with an invariant mass 
within 25 MeV/c 2 of the K° s mass as given by the Particle Data Group (PDG) 497.672 MeV/c 2 [12] . 
Only events with exactly one Kg candidate are retained. The track from the signal side not 
originating from the Kg candidate is required to be identified as a pion and originate from the 
interaction point. Pions are identified by the ionisation loss in the tracking system (dE/dx), the 
shape of the shower in the calorimeter and their discrimination from kaons performed in the DIRC. 
All tracks on the signal side are required to lie within the geometrical acceptance region of the 
EMC and DIRC to ensure good particle identification. 

Additional cuts are imposed to further reduce the backgrounds. The net charge of the event is 
required to be zero and a cut requiring the thrust of the event to be greater than 0.85 is imposed 
to reduce the non-r background. 

Backgrounds from Bhabha events are suppressed by requiring the momentum of the lepton-side 
track to be less than 4.9 GeV/c. Backgrounds from radiative Bhabha and //-pair events with a 
converted photon are suppressed by requiring the modulus of the cosine of the decay angle to be less 
than 0.97. The decay angle is defined as the angle between the momentum of the ir + originating 
from the Kg in the Kg's rest frame and the Kg momentum in the laboratory frame. When this 
quantity is calculated on electrons misidentified as pions the effect of assigning an incorrect mass 
will push the value towards ±1. From studies of missing transverse event energy, backgrounds from 
two-photon events are determined to be negligible. 

Along with signal events, r~ — > ir~ir + ir~v T events also produce three pions in the final state, 
but all three come directly from the primary interaction point. To remove r~ —* ■K~ir + ir~i> T events 
the Kg flight length significance in the plane perpendicular to the collision axis is required to be 
greater than 5.0. The flight length significance is defined as the measured flight length divided by 
its estimated uncertainty. This cut removes approximately 90% of the remaining r _ — > tt^tt + tt^u t 
events in the sample. FigureQ] shows the Kg flight length significance distribution for events in this 
analysis with all other cuts applied. 

To increase the likelihood that the pions from the Kg candidate really come from a Kg they 
are required to have a distance of closest approach to each other of less than 2 mm. The Kg 
trajectory is reconstructed by performing a vertex fit on the Kg daughter pions constraining them 
to originate from a single point and then summing their momenta. To increase the probability that 
the Kg originated from the interaction point this trajectory is required to have a distance of closest 
approach from the collision axis of less than 1 mm. 

Once these cuts have been applied, the largest two background channels are from t~ — > 
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Flight length significance L /a 



Figure 1: The K% flight length significance in the zy-plane for the the combined (e-tag+/x-tag) sample. The MC 
signal branching fraction is set to the measured value from this analysis. 



K^KgTT~ v T and t~ — * Kgir~ir °u T events where the additional neutral particle to the Kg is unde- 
tected. In order to reject such events maximum neutral energy cuts are applied. The total energy 
in the calorimeter not associated to any charged track must be less than 0.5 GeV. In the signal 
hemisphere this quantity is required to be below 0.25 GeV. 

After selection about 80% of the retained MC events are r — > K s tt~u t . About 98.5% of the 
background events come from r-decays and about 1.5% from non-r-events. The overall signal 
efficiency is estimated from the MC simulation as: 

/yscl 

e ■ = 3ig (1) 

sig 

where AJ™ is the number of generated signal events and A^ 1 is the number of selected signal events. 
A total of 5.872 x 10 6 r — ► A°7r~z/ T signal events were generated. 

The total selection efficiency, which includes the efficiency corrections that are described in 
Section HJ is (0.620 ± 0.003)%, (0.482 ± 0.003)% and (1.101 ± 0.004)% for the e-tag, //-tag and 
combined samples respectively. 



4 EFFICIENCY CORRECTIONS 



Since imperfect detector simulation may mean that the reconstruction/selection efficiencies differ 
between data and Monte Carlo, some corrections are applied to Monte Carlo signal and background 
events. 
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A detailed study of the K s efficiency has been carried out. B° — ► 7r + D~ (D~ — > K s tt~) decays 
were used to study the K s efficiency for data and Monte Carlo. As a result of this study, a per K s 
average data/MC efficiency correction of 0.983 should be applied to the MC sample. Since only one 
K^ is selected in this analysis the MC sample is weighted with 0.983 as a Kg efficiency correction. 

The efficiency of particle identification (PID) cuts are known to differ between data and Monte 
Carlo and so corrections to the Monte Carlo are applied. A set of efficiency tables, binned in p, 9 
and (p, are used in order to obtain the necessary weights to use to correct the Monte Carlo events. 
The efficiency correction (data/MC relative efficiency) and corresponding uncertainty is calculated 
for each track. 

In this analysis charged PID selection is performed on the pion not originating from the Kg 
and the lepton. The average values obtained for the corrections due to the electron, muon and pion 
are 1.015, 0.948 and 0.979 respectively. 

The total data/MC efficiency correction, e corr , is made by combining the efficiency corrections 
described above. e corr is used to weight the MC sample and the average values obtained are 0.977, 
0.912 and 0.947 for the e-tag, /z-tag and combined samples. The systematic error coming from this 
procedure is described in Section [5j 

5 SYSTEMATIC UNCERTAINTIES 

A number of systematic uncertainties have been considered. These are described below. 

An error of 1.10% is included due to the K s efficiency correction procedure described in Sec- 
tion HI When propagated through to the branching fraction measurement this gives a 1.40% 
systematic uncertainty on the branching fraction. 

The tracking efficiency is susceptible to bias caused by physics data and MC simulation dis- 
crepancies. An error of 0.23% per track not originating from the Kg is assigned to account for 
this. Tracking uncertainties from the K s daughters are included in the K s selection systematic. 
Each event contains two reconstructed charged tracks with correlated uncertainties, leading to a 
total tracking efficiency uncertainty of 0.46%. This gives a 0.58% systematic uncertainty on the 
branching fraction. 

The total particle identification uncertainty that arises from the efficiency correction procedure 
described in Section H] for the e-tag, //-tag and combined samples is estimated to be 1.05%, 1.33% 
and 1.18%, respectively. When propagated through to the branching fraction measurement this 
gives a systematic uncertainty of 1.45% for the e-tag, 1.68% for the //-tag and 1.50% for the 
combined sample. 

The relative uncertainty associated with the t + t~ pair production cross-section and the BABAR 
luminosity determination is 0.65%, leading to a 0.83% systematic uncertainty on the branching 
fraction. A relative error is applied for uncertainty in the modelling of selection variables on which 
cuts are applied. The modelling efficiency uncertainty is estimated as 0.29%, leading to a 0.37% 
systematic uncertainty on the branching fraction. 

The systematic uncertainty on the signal efficiency and branching fraction measurement due to 
signal Monte Carlo statistics is 0.51% for the e-tag, 0.66% for the //-tag and 0.37% for the combined 
samples. The error on the number of background events due to limited Monte Carlo statistics is 
1.02% for e-tag, 1.12% for //-tag and 0.76% for the combined sample. This gives a systematic 
uncertainty on the branching fraction measurement of 0.28% for e-tag, 0.30% for /i-tag and 0.20% 
for the combined sample. 
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Since a number of r decay modes have not yet been precisely measured, particularly those 
modes which contain Cabibbo suppression factors, the branching fractions used as input to the 
Monte Carlo simulation come with an uncertainty which feeds into the total systematic error. In 
order to evaluate this a weighted-sum of the r backgrounds is constructed using the Monte Carlo 
truth for r Monte Carlo events passing the analysis selection criteria. 

The branching fraction for the background mode t~ — ► K~ K + tt~u t is taken from a recent 
BABAR collaboration measurement B(t~ K~K + tt~u t ) = (0.1346 ± 0.0010 ± 0.0036)% [14J. We 
predict the ratio £>(t~ — * tt~ 'KgK^i/ T ) /B(r~ — > K~K + tt~v t ) to be 0.50 ± 0.05 using an isospin 
relation [17j . This gives B(r~ -> tt-K%K\v t ) = (0.0673 ± 0.0070)% which is consistent with the 
PDG 2006 value B[r~ -> -k~K%K\v t ) = (0.112±0.030)%. All other background mode uncertainties 
are taken from the PDG 2006 [12]. The overall estimated uncertainty due to the t backgrounds is 
given by: 




where Wi is the fraction of selected background r events of mode i, B t is the branching fraction of 

mode % and a t is the uncertainty of Bi. 

Table[T] shows the weights and uncertainties of the r background modes remaining in the selected 
— ► Kgir^Uj- sample. The resulting uncertainty attributed to the r background modes (A T_bkg ) 

is estimated as 4.99% on the number of background events and is consistent between each of the 

tagged samples. This leads to a 1.37% systematic uncertainty on the branching fraction. 



Table 1: r background uncertainties and their weights in the selected r — > K%-k v t sample. 



Decay Channel 


w[%] 


m 


t~ -> ir-K^v T 


24.24 


10.53 


t~ -> K-K° s TT a u T 


0.57 


17.54 


t~ tt-K°K° l v t 


41.16 


10.38 


t~ -» ■k~K°K q s v t 


0.18 


20.83 


T~ -»■ K-R°u T 


2.64 


10.46 


T~ — > ir~Tr°v T 


0.75 


0.39 


T~ — > 2h + tt~v t 


23.65 


0.89 


T~ — y 27T~ 1T + 7T° V T 


3.41 


1.35 


t~ —> K~ K + ir~i> T 


0.61 


2.78 


T~ — > K~ir + ir~ v r 


0.78 


10.26 



Table [2] summarises the main sources of systematic uncertainty in the analysis. 

The relative systematic uncertainty on the branching fraction includes the individual sources 
described above when propagated through equation ([3]) taking into account correlations between 
the different components considered. The total systematic uncertainty on the branching fraction 
measurement is 2.73%, 2.87% and 2.72% for the e-tag, //-tag and combined samples, respectively. 

6 BRANCHING FRACTION MEASUREMENT: £(r K°tt ~v t ) 

Figure[2] shows the invariant mass spectrum of the selected data and MC K^tt~ candidates in the 
combined (e-tag+/i-tag) sample after all the analysis requirements, including efficiency corrections. 
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Table 2: Summary of the systematic uncertainties as they feed into the measurement of £>(r — > K°ir v T ). 



Systematic 


e-tag 


//-tag 


Combined 


Tracking 


0.58% 


0.58% 


0.58% 


K% Efficiency 


1.40% 


1.40% 


1.40% 


PID 


1.45% 


1.68% 


1.50% 


C x a TT 


0.83% 


0.83% 


0.83% 


Statistical efficiency error 


0.51% 


0.56% 


0.38% 


MC background statistics 


0.28% 


0.30% 


0.20% 


r backgrounds 


1.37% 


1.37% 


1.37% 


Modelling efficiency 


0.37% 


0.37% 


0.37% 


Total 


2.73% 


2.87% 


2.72% 



In this plot, the signal MC sample is scaled using the branching fraction measured in this anal- 
ysis: B(t~ — ► Kgir^Ur) = 0.420%. The data and MC sample disagreements can be ascribed to 
inadequate modelling of the K* resonances in the MC simulation. There is some evidence that the 
K*(892) mass has been underestimated in the MC simulation and of the existence of Kq (800) and 
Ko(1430) resonances which are not modelled [7]. However these discrepancies do not affect our 
result for the branching ratio. 

The final result for B{t~ — > K®ir~v T ) is estimated using the combined e-tag + /^-tag sample, 
where the total number of events observed, estimated background level and efficiency are derived 
from the two tagged samples. As a cross-check, B{r~ — > K°tt~v t ) is calculated for each tagged 
sample separately and the results obtained are in excellent agreement with the combined result. 

The branching fraction B(t~ — > K®n~v T ) is estimated by 

B{r- - KH-u T ) = 1 N ^~ N ^^ (3) 

TT ^Sig 

where N TT is the total number of r + r _ pairs in the real data, -/V data is the number of selected events 
in real data, -/V bkg is the number of background events estimated from Monte Carlo sample and e sig 
is the signal efficiency as defined in equation ([T|). The total number of t + t~ pairs in the data is 
given by 

N TT = <7 T £ data = (353.4 ± 2.3) x 10 6 , (4) 

where a T is the t + t~ production cross-section at BABAR (i.e. 0.919 ± 0.003 nb) and £ d at a is the 
(integrated) real data luminosity (i.e. 384.6 ± 2.2 fb -1 ). 

Table [3] gives the numbers of real and simulated data events passing the selection criteria that 
feed into the £>(r~ — > K°tt^u t ) calculation. A result for the measurement of the branching fraction 
B(t~ — > K q tt~v t ) is given in Table HI 

Figure [3] shows previously published measurements of B(r~ — > K Q/ n~v T ). It can be seen that 
this BABAR preliminary result is the world's most precise measurement. 
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Figure 2: Reconstructed r~ — > iC'gTr - ^ mass distribution for the combined (e-tag+/x-tag) sample using (left) 
logarithmic scale and (right) linear scale. The MC signal branching fraction is set to the measured value from this 
analysis. 
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Figure 3: Previously published measurements of £>(r — » K°ir v T ). The Belle result and the result from this 
analysis do not contribute to the PDG 2006 average. The BABAR preliminary result provides the world's most precise 
measurement to date: £(t~ -> K°Tr~v r ) = (0.840 ± 0.004 (stat) ± 0.023 (syst)) %. 
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Table 3: Numbers of data and Monte Carlo events remaining in the selected samples corresponding to the 384.6 fb 
(integrated) real data luminosity. The signal Monte Carlo has been scaled with the PDG 2006 value |12j . 



Data e-tag A*-tag Combined 

Real" 47092 ± 217 36641 ± 191 83733 ± 289 

Signal MC 39445 ± 193 30749 ± 176 70194 ± 261 

t + t- background 9942 ± 92 7645 ± 88 17587 ± 131 

uds 8.9 ± 3.4 65.2 ± 7.9 74.1 ± 8.5 

cc 45.5 ± 5.0 43.5 ± 4.8 89.1 ± 6.9 

BB 3.4 ± 1.1 3.7 ± 1.2 7.1 ± 1.6 

fjt+fj,- ± 14.2 ± 3.7 14.2 ± 3.7 



Table 4: BABAR preliminary 6(r~ -> K" 7v v T ) measurements. 
Sample g(r~ -> K^~v T ) [%} 
e-tag 0.840 ± 0.005 (stat) ± 0.023 (syst) 

//-tag 0.840 ± 0.006 (stat) ± 0.024 (syst) 

Combined 0.840 ± 0.004 (stat) ± 0.023 (syst) 



7 SUMMARY 

Using 384.6 fb -1 of e + e~ collision data produced by the PEP-II collider and recorded by the BABAR 
detector, we obtain the preliminary result: 

(0.840 ±0.004 (stat) ±0.023 (syst) ) %. (5) 

This result is the most precise measurement of this branching fraction to date and is consistent 
with the world average. 
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